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Abstract 

The process of formation of fractal structure in one-dimensional 
self-gravitating system is examined numerically. It is clarified that 
structures created in small spatial scale grow up to larger scale through 
clustering of clusters, and form power-law correlation. 



1 Introduction 

l-dimensional self-gravitating system has long been used to study systems 
interacting through self-gravity [Por32| , ^lF67| , |RM88| , |YGS92| , Ircfel . The 



model helped us to understand various properties, such as motion of starts 
normal to the galactic plane, relaxation properties of elliptical galaxies, and 
so on. 

Recently we have found that fractal structure dynamically emerges in the 
system | KKUCI|| . It is surprising that the fractal structure is created from 



non-fractal initial conditions. 
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The emergence of fractal structure has several important meanings. First, 
there are many fractal (or power-law type) structure in nature. In astro- 
physics the distribution of young stars, galaxies and cluster of galaxies are 
known to obey power- law [ [Lar95| , |NTHN98| , PCGE"^9^ . However there is 
no universal explanation about the origin of power-law behavior in general. 
Since the dynamical emergence of fractal structure arises in self-gravitational 
model, there can be a relation between the origin of these fractal structures. 

Second. There are many studies about the relaxation process of 1-dimensional 
self-gravitating system and it is reported that relaxation is quite slow ( ||RM88| , 
l'GK94]|, and references therein). It is interesting if there is any relation be- 



tween the slow relaxation in the model and the dynamical emergence of 
fractal structure. 

Since the system is Hamiltonian system with many particles, one may 
expect that most initial condition lead the system to thermal equilibrium. It 
is interesting that some initial conditions, e.g., uniformly random in spatial 
distribution, evolves to form fractal spatial structure, instead of just become 
thermalized. Hence to study the origin of the dynamical emergence of fractal 
structure may give unique insight to non-thermal behavior and structure 
formation in many-body systems in general ||KK92| , |KK94 ] . 

Therefore it is quite interesting to investigate the reason why the fractal 
structure emerges from non-fractal initial conditions in 1-dimensional self- 
gravitating system. In order to do that, it is essentially important to know 
in detail the process how the fractal structure is created, not just to observe 
the grown-up fractal structure. 

In this letter we continue the study to investigate the process how the 
fractal structure is developed through time evolution. 

The model is presented in section 2. In section 3 we will examine in detail 
how the power-law correlation is formed, and the final section is for summary 
and discussions. 



2 Model 

The model we use in this letter is the same as the one used in our previ- 
ous letter ||KKOO|| , that is, the one- dimensional self-gravitating system. The 
Hamiltonian is 



N 



H 



i=l 



+ \xi — xA , — oo < Xj < oo . 

2m ^ 



i>j 



This model represents one-dimensional motion of parallel sheets of infinite 
extent interacting through Newtonian gravity. 
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The model is simple and tractable, and yet it also contains essential fea- 
ture of gravitational interaction. 
Throughout this paper we set 

m = 1/N and AttG = 1 . (2) 



3 Formation of power-law correlation 



In the previous letter |[KKOO|| we showed that fractal structure emerges from 
non-fractal initial conditions. Typical initial conditons we used are those 
of virial ratio 2Ekin/Epot = 0. (Spatial distribution is set to be random.) 
These state of zero velocity dispersion corresponds to the limiting case of 
zero thermal fluctuation. 

To analyze the spatial structure we use two-point correlation function 
^{r), which is defined as 

dP = ndV{l + ^{r)) (3) 

where dP stands for probability to find another particle in volume dV at 
distance r from a particle, n is the average number density. If particles are 
distributed independently and no spatial structure is formed, then dP does 
not depend on r and ^{r) = 0. If, on the other hand, ^(r) increases as r 
decreases, the probability that we find another particle near one particle gets 
large, and we can understand cluster-like structures are formed. If the ^(r) 
does not have characteristic length, then we can understand that the cluster 
system does not have characteristic length scale. 

Fig.0 shows temporal evolution of ^(r) for early stage of time. We can 
clearly see how the power-law behavior of correlation function grows up. 
That is, power-law structure is first created in small scale, then gradually 
grows up to larger scale. 

The graphs are aligned from bottom to top as time proceeds because 
we are looking at small-sized scale here (r < 0.002 ) and particles initially 
separated are coming near, which increase ^(r) for small r. 

The time dependence of two point correlation function C, {r) tells us when 
and where the power-law structure is created and how it grows up. However 
C,{r) does not tell what kind of structure is actually formed. So we need 
another quantity to understand the detail of the structure formation. 

Since we are dealing with a gravitational system, interaction is attractive. 
Each particle, initially distributed in the 1-dimensional space, attracts each 
other by gravitational force, and small clusters are formed here and there. 
These clusters attracts and moves around each other, as well as their internal 
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two-point correlation function xi(r) 




Figure 1: Evolution of correlation function ^(r) at t = , £ = 3, 4, ■ • ■ , 17 
(from bottom to top). System size N = 2^^, Initial condition: Xi =random, 
Vi = 0. Etotai = 0.25. 



oscillation (rotation, in 3D). While these processes go on, the order of parti- 
cles given at initial condition is exchanged, and initial distribution, aligned 
on the X-axis as seen in the {x,v) phase space (/i-space), get folded many- 
times. Hence the degree of folding is a good quantity to understand how the 
cluster structure is formed. 

We can measure the number of folding as follows. Since we are considering 
the case where initial conditions have zero velocity dispersion, distribution 
of particles in {x, v) space (/i-space) in each initial condition is on a curve, 
and we can label particles along the curve. For example, if we take the initial 
conditon as 

Vi{t = 0) = 0, « = 1,2,--- ,iV , 
we can label the particles as 

xi(0) < X2(0) < ■■■ < xn{0) . (4) 
Then we can count f{t), the number of folding at time t, as 

f{t) = number of i {1 < i < N) which satisfies 

{xi+i{t) - Xi{t)) ■ {xi{t) - Xi_i{t)) < (5) 

Fig. 1^ shows how the initial distribution becomes folded. Time depen- 
dence of number of folding f{t) shown in this figure and time dependence 
of the two point correlation function ^ (r) shown in Fig. |l] corresponds well. 
That is, power-law behavior of ^(r) develops from small scale as the number 
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Figure 2: Time evolution of number of foldings. Horizontal axis represents 
time measured in steps i defined in eq.(|) 



of folding increases, and become most clear when f{t) is at a maximum at 

t = tlQ. 

The process of structure formation, as seen in two quantities ^(r) and 
f(t), can also be understood in the real space of particle distribution. Fig. |^, 
represents time evolution of particle distribution of particles {xi{t),Vi{t)) , 
8200 < i < 8300 for t = ts,ti3, tis, and tiQ, where ti are time steps defined as 

- |i< • (6) 

(These time steps are defined just for practical purpose.) 

In these figures we connect each particles on {x, v) space according to the 
order defined at the initial condition (g) so as to enhance the structure of 
folding. Vertical steps on the graphs represent the location where particles 
are clustered. 

At t = ts we can see several small clusters are formed. These clusters 
gradually approach each other to form larger clusters, as seen in t = tis, tis 
and ti6 in the figure. Through this process power-law correlaion is developed. 

In later time the clusters get together and larger clusters are developed. 
In the left of Fig. ^ we show a snapshot of the distribution of particles shown 
in the previous figure Fig.|[ This cluster itself is a part of a larger cluster 
system shown in the right of the figure. 
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Figure 3: Distribution in {x,v) space at t = tg (up:left), ti3(up:right), 
ii5(down:left) and tie (downiright) for particles No. 8201 ~ 8299 . Par- 
ticles are linked according to the order assigned at t = 0. 
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Figure 4: continued from Fig.^ (left) and /i-space (wide ranged) (right) at 

t = t55 
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4 Summary and discussions 



In this letter we clarified the process of formation of power-law correlation in 
one-dimensional self-gravitating system, found in our previous letter |[KKOO| . 
Power-law behavior of two-point correlation function begins in small spatial 
scale, then grows up to larger scale. Number of foldings in the order of parti- 
cles gets large when the power-law behavior is clearly seen. These results and 
phase space portraits represents the process of formation of fractal structure 
as follows: tiny clusters are formed first, then the clusters get together to 
form larger clusters, and so on. 

The process we showed is quite similar to the "hierarchical clustering 
PeeSOj , [Man77|| , which is used for descibe the large-scale structure 



scenario 



of universe. It is known that two-point correlation function of galaxies and 
cluster of galaxies obey power-law |pCGE'^"94[ . An explanation for the origin 



of the power-law structure is that density fluctuation of small scale is created 
first to become some structure of small scale (bottom- up), and, according 
to the "hierarchical clustering scenario" the structure made in small scale 
attract each other through gravitational interaction and merge into larger 
structure, and so on. The initial state is now considered to be of small 
temperature fluctuation. 

In this letter we have shown one example which obeys this "hierarchical 
clustering scenario" by purely gravitational interaction. An important fact is 
that our result shows that power-law correlation emerges through hierarchical 
clustering. Also the initial condition we found for fractal structure to be 
created is the zero velocity dispersion, which coinsides with zero temperature 
fluctuation considered in the cosmology. 

The clustering process is developed by gravitational interaction. Hence 
similar process can be actually occuring in nature where the dynamics of the 
system is governed by gravity. One candidate is the power-law distribution 



of yound stars in star-forming regions [Lar95, NTHN98 



The emergence of power-law correlation can be considered as a subject 
of structure formation in general. Our result shows a process of generating 
fractal structure by purely dynamical interaction. 

In this letter we have clarified how the fractal structure is formed from 
non-fractal initial conditions in one-dimensional self-gravitating system. We 
hope that this result helps us to understand why the structure is created. 
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